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Abstract 

We study the Von Neumann and Renyi entanglement entropy of 
long-range harmonic oscillators (LRHO) by both theoretical and nu- 
merical means. We show that the entanglement entropy in massless 
harmonic oscillators increases logarithmically with the sub-system size 
as S — %p-logL Although the entanglement entropy of LRHO's 
shares some similarities with the entanglement entropy at conformal 
critical points we show that the Renyi entanglement entropy presents 
some deviations from the expected conformal behavior. In the massive 
case we demonstrate that the behavior of the entanglement entropy 
with respect to the correlation length is also logarithmic as the short 
range case. 



1 Introduction 

While the entanglement entropy of short-range interacting quantum systems 
in one dimension is well studied with many different techniques pQ, such as 
hamiltonian techniques [2H5], euclidean methods [HIE] and conformal field 
theory [8l[9] there are few results concerning the long-range interacting sys- 
tems. The main difficulty is the lack of exact solution for these systems which 
makes the problem much harder than the short-range counterparts. Because 
of the presence of non-trivial dynamical exponents the common euclidean 
techniques are usually useless in calculating the entanglement entropy in 
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these systems. It is worth mentioning that although having a non-trivial 
dynamical exponents is not restricted to just long-range interacting systems, 
see for example [lOy ilj , they usually provide a very natural knob to change 
it arbitrarily. Because of the huge finite size effects in these systems even 
the numerical calculations are not easy as their short-range counterparts. 
Nevertheless, recently much progress has been made in different directions: 
In [12j the entanglement entropy in Lipkin-Meshkov-Glick (LMG) model is 
studied. This model has hamiltonian similar to the XY model but while in 
the latter model the interaction only takes place between nearest neighbors, 
in the LMG model, all spins interact among themselves. In [13], the interac- 
tions are restricted to the Ising-type without external magnetic field which 
allows to study both the static and the dynamical entanglement properties 
of the system. Eisert et. all [14] found a logarithmically divergent geomet- 
ric entropy in free fermions with long-range unshielded Coulomb interaction. 
Plenio et.all [15] studied the general properties of the entanglement entropy 
for interacting harmonic oscillators. In particular they found that the area 
law should be valid in higher dimensions even for long-range interacting har- 
monic oscillators, and most recently the entanglement entropy for the long 
range anti- ferromagnetic Ising chain is calculated by numerical means [16] . 

Harmonic oscillators are the building blocks of many quantum mechan- 
ical and field theoretical systems and of course this is also true for most of 
the long-range interacting quantum systems and non-local field theories. In 
this work we study the entanglement entropy of long-range harmonic oscil- 
lators in the massless and massive cases. The main approach we have used 
is the so called hamiltonian technique, which was first introduced in [2] and 
then elaborated in [7] , for review see [5j[6] . The important advantage of this 
method with respect to the others in calculating the entanglement entropy 
in our problem is that, firstly it is possible, to a large extent, to carry out the 
calculations exactly, and on top of that the numerical calculations are very 
easy to implement. We will consider the lattice and continuum calculations 
in parallel. We first introduce the method and then apply it directly to our 
problem. In the analytical side we will write the main eigenvalue problem ex- 
actly, then we will try to find approximate solutions by using our numerical 
calculations. Using the solutions of the eigenvalue problem we then also cal- 
culate the Renyi entanglement entropy and finally we will solve numerically 
the problem of the massive interacting long-range harmonic oscillators. 
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2 Results 



The entanglement entropy of a subsystem A is denned by using reduced 
density matrix p A as 

S A = -tip A logp A - (1) 

In this work we consider three main kinds of subsystems. In the massless 
case we consider two kinds of systems: in the first one the system is large 
and A is a small sub-system with length I and in the second case the system 
is finite with length L and the sub-system is just half of it with length L/2. 
In the massive case the system is large and the sub-system is just half of 
the system. We will provide approximate formula just for the first case and 
study the two other cases numerically. 

We define the hamiltonian of one dimensional interacting harmonic os- 
cillators as 

TV N 



2 ^ T 2 

i=i jj=i 



where on the lattice we take 



m 



r(-f + i-j)T{a + l) . a a 
irlXl + S + i-j) s ^r ) + m ' 



(3) 



In the special case when a/2 is an integer, K(n) = (— l) a ~ n+1 C a ^ +n , 
where C a ^ +n are binomial coefficients. In this case we remark that K(n) = 
for n > a/2. It is easy to see that in the special case a = 2 the K matrix 
is just a simple laplacian. For non-integer values for large distances we have 
K(n) ~ 

In the continuum limit the above hamiltonian can be written as 

TV 

\ £ hKirfj -> f [-^(x)(- V ) Q/2 0(x) + l -m^]dx, (4) 

where — (— v) a ^ 2 is the fractional laplacian defined by its Fourier transform 
\q\ a . 

To measure the entanglement entropy of a sub-system with length I of 
an infinite system we follow the method explained in [7]. First we make the 
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matrices W = K 1 / 2 and W 1 = K 1 / 2 . In the continuum limit they have 
the following forms 

W ±x (x,y) = — f {\k\ a + m a ) ±l l 2 e ik -( x - y) 

11 

= l C(m a ) 

2T[Tu/2]cos(^)(x-y) 1±a / 2 1 ; ' 
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Figure 1: (colour online). Entanglement entropy for ID LRHO: entan- 
glement entropy changes linearly in the semilog plot in the region < I < 
L/100, here we have L = 4000. From top to bottom a = 1.8, 1.6, 1.4, 1.2, 1.0. 

Then we define the matrix A by multiplying W and W~ l in the sub- 
region O £ (— oo < z < 0)U(l < z < oo) that we wish to calculate its entan- 
glement entropy 



A(x,y) 



A 



dz[W- 1 {x,z)W(z,y) 



l—x 

i-v 



a/2 



a/2 



a(x - y) 



(6) 



where A - 
ir 2 A(x,y) 



4T[-a/2]T[a/2] cos 2 (^)- For Q 



2 the matrix has a different form 



{1 - X)l %^:^ 0g[l - V) - Ei2S gEgjr SM - We campaired the 
matrix A(x, y) coming from the above equation with numerical A and found 
a very good agreement when the distances are more than four lattice sizes. 
The agreement gets better by increasing the size of the system. 
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The entanglement entropy can be expressed directly in terms of the 
eigenvalues E{ of the matrix A as [21[7] 

S = 5>g ^ + v^T^log(-^ + + (7) 

Ei 

In the continuum limit one needs to solve the following eigenvalue problem 

J ' dyA(x,y)1>(y)=Eil>(x). (8) 

At the numerical level we followed the above algorithm and calculated 
the entanglement entropy for large enough system sizes such as L = 5000 and 
L = 6000, to avoid any finite size effects [T7]. For these sizes the data show 
stability even for very small a's as far as we take the sub-system size less 
than L/100. Interestingly the entanglement entropy grows logarithmically 
with the sub-system size as 

S = ^logl, (9) 




Figure 2: (colour online). Green (Dark gray) circles represent c e //(a) for 
discrete systems with size L = 6000. The prefactor c e ff(a) is measured 
using the scaling relation S with log I by Chi-square test analysis in the 
region < I < L/100. The red line represents the same quantity coming 
from the continuum limit approximation. The prefactor c e //(a) = 3j^| 
increases monotonically with increasing a and it is equal to the central 
charge of the free field theory at a = 2. 
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where c e ff is a dependent, see Fig 1. The c e ff is maximum at the 
conformal short-range point and it starts to decrease by decreasing a. For 
a = 2 Callan and Wilczek [7] found a good estimation for the eigenvalues 
of the A operator when the subsystem is half of a finite system. Motivated 
by their work we found E{uj) = sinh ^ nuJ ^ , where oj(Ei)\og(l) = ^ a very 
good approximation for the eigenvalues of the short-range case. Based on 
numerical comparison we found that apart from a constant the behavior of 
the logarithm of the small eigenvalues (large i, i.e. i > 6) is independent of 
a and one can safely conjecture the following behavior for the eigenvalues 
of the A operator 

E(oS) = , (10) 

sinh (ttu) + 6(a) 

where a(a) and 6(a) are functions of a and wlog(Z) = ^. Non-zero value 
of 6(a) is necessary to make the formula compatible with larger eigenvalues. 
We found a{a) = § sin 2 (^) and b(a) = 0.1146a + 0.1868a 2 - 0.1988a 3 + 
0.03771a 4 best fit functions to our numerical data. The value of 6(a) is 
zero at a = and 2 and it has a maximum at a = 1. In principle 6(a) 
might be not a simple polynomial. Using the above formula one can get the 
logarithmic behavior for the entanglement entropy for free and the c e fj is 

ceff(a) = - [°° du\\og^^ + ^l + E{u) 

n Jo L z , , 

1 / — - 

X log( == + Jl + 



In Fig 2 we compare the result coming from the above formula and the 
numerical calculations, with excellent agreement. The logarithmic behavior 
of the entanglement entropy is reminiscent of the presence of conformal sym- 
metry. To check this point we calculated the entanglement Renyi entropy 
using the same technique. The Renyi entropy S n is defined as 

S n = -^—logtrp n , n>\ (12) 
1 — n 

and can be calculated using the folowing formula [6] 

S n = V (log(l - e-^) - nlog(l - e^)) . (13) 

n — 1 — ' 



where e £i = 1 + — \f\ + El) with Ei as the eigenvalues of the A matrix. 
In conformal invariant systems S changes logarithmically with a prefactor 
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Figure 3: (colour online), versus n for different a's (from top to bottom 
a = 2.0, 1.8, 1.6, 1.4.). The prefactor c^j(a) is measured using the scaling 
relation of S n with log/ by Chi-square test analysis in the region < I < 
L/100 (The system size L = 6000). Red solid lines come from the continuum 
limit approximation using the equation ()10p . Inset: A and B coefficients 
versus a for discrete systems with size L = 6000. At the conformal point 
they are equal to one. 



|(1 + -), where c is the central charge [9]. We found that in LRHO the 
Renyi entropy is also increases with the sub-system size logarithmically as 

S n = -^-log(Z). We have depicted c^j-j versus n in Fig 3. The best fit is 
c™jj = ^Tp-(^4(a) + B(a)/n) where A(a) and B(a) are functions of a. The 
a dependency of these coefficients shows that this system is not conformally 
invariant except at a = 2 where we have A = B = 1. 
In the continuum limit one can write (|13|) as 

Sn = -^\-, /°°^(log(l-C)-nlog(l-0), (14) 
vr(n - 1) Jo 



where £ = v . — and E(uj) comes from Eq. (flUI) , After doing the 

integration numerically one can find c"^j versus n, see Fig 3. The perfect 
agreement between numerical calculation and the continuum limit results is 
another support for the accuracy of the equation (jlOp . 

Now we discuss the effect of different boundary conditions to the entan- 
glement entropy. The most interesting one is the free boundary condition, 
which we take a finite system with half of it as the subsystem. In this case 
we define the K matrix by throwing away the elements of the infinite ma- 
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Figure 4: (colour online). The blue (black) stars are the prefactors of the 

logarithm for the massless system with length L, i.e. -^p- = <S/log.L. The 
system size L = 1000 — 4000 used to measure crf^ by the best fit to S 
versus logL. The light green (light gray) circles are the prefactors for the 
entanglement entropy of the massive systems with the system size L = 4000, 

c 9 

i.e. = S/ log m. The quantity c g e ^ is measured using the best fit to 5 in 
the saturation regime versus logm by chi-square test (m G (lO -3 — 10 -2 )). 

trix which is not inside the system. The immediate problem that we face 
is then that we are not able to write the W matrices analytically because 
in the presence of the boundary we are not allowed to use Fourier trans- 
form. In this case we follow all the calculations numerically. The algorithm 
is the same and interestingly we found again logarithmic behaviour for the 

entanglement entropy with respect to the size of the system -^p logL. The 
prefactor c^f* is in general different from c e // except at conformal point, 
see Fig 4. This fact indicates that in the long-range interacting systems the 
entanglement entropy of a sub-system is not equal to the summation of the 
contribution of each boundary of the sub-system with its complement. 

Finally we report the entanglement entropy of massive long-range inter- 
acting harmonic oscillators. The entanglement entropy of the positive half 
of an infinite massive system for short-range models comes from the Cardy- 
Calabrese formula S = | log m where c is the central charge of the system 
and it is one for short-range harmonic oscillators [9] . We calculated the same 
quantity for the long-range interacting systems. We found that S saturates 
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in I — > oo limit, and changes logarithmically with respect to the mass as 

c 9 

S=^logm, (15) 

where cP e ^ is different from c e ff except at the a = 2, however, interestingly 
it is equal to the prefactor of the free boundary condition case, see Fig 4. 
The nature of this equality is not clear for us. Since in long-range interacting 
systems the correlation length £ changes as m a one might expect that the 
the entanglement entropy be proportional to log£ and then expect that 
c e// ^ e proportional to a. Although it seems that this is the case in the 
region 1 < a < 2 the slop that we find is not compatible with the natural 
expectation c? e ^ = at a = 0. Moreover our primary numerical calculations 
show strong deviation from the linear behaviour close to a = 

3 Conclusion and outline 

: In this paper we have found the entanglement entropy of long-range har- 
monic oscillators using the hamiltonian technique. When we consider a 
sub-system with length I the entanglement entropy follows the logarithmic 
behaviour as the short range cases but the prefactor is dependent on the 
power of the interaction a. For stronger long-range interactions the pref- 
actor is smaller. Although it is tempting to interpret the prefactor of the 
logarithm c e ft as the effective central charge of the system, by calculating 
the Renyi entropy we showed that the nature of this prefactor is differ- 
ent from the central charge at least at the level of the Renyi entropy. We 
also calculated entanglement entropy for two other interesting cases, the 
system with boundary and also the massive case. Interestingly we found 
that the prefactor is the same in these two different cases, however, it is 
different from the massless infinite system. Our study shows that the entan- 
glement entropy of long-range interacting systems shows some similarities 
with the entanglement entropy of short-range systems , however, since they 
have very different symmetries they start to show strong differences when 
we study their replica behaviour. We believe that long-range interacting 
systems show interesting features which deserve more intense studies. Our 
work can be extended in many directions: calculating the entanglement en- 
tropy in higher dimensions, the dynamics of the entanglement entropy and 
the finite size effects are just few among many other interesting directions. 
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